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DETERMINATION OF MEAN CAMBER SURFACES FOR WINGS HAVING UNIFORM 
CHORDWISE LOADING AND ARBITRARY SPANWISE LOADING IN SUBSONIC FLOW^ 

By S. EIatzoff, M. Fhancbs Faison, and Hugh C. DuBosb 


SUMMARY 

The field of a uniformly loaded wing in subsonic flow is dis- 
cussed in terms of the accderalion ^potential. It is shown thal, 
for the design of such wings ^ the slope of the mean camber surface 
at any point can be determined by a line integration around the 
wing boundary. By an additional line integration around the 
wing boundary f this method is extended to include the case where 
the local section lift coefiicient varies with spanwise location {the 
chordwise loading at every section still remaining uniform). 

For the uniformly loaded wing of polygonal plan form, the 
integrations necessary to determine the local slope of the surface 
and the further integration of the slopes to determine the ordinate 
can be done anxdytically. An outline of these integrations amd 
the resulting formulas are mduded. 

Calculated results are given for a sweptback wing with uniform 
chordwise loading and a highly tapered spanwise loading, a 
uniformly loaded delta wing, a uniformly loaded sweptback wing, 
and the same sweptback wing with uniform chordwise loading 
but elliptical span load distributio'n. 

INTRODUCTION 

The design of mean camber surfaces to sustain a specified 
area distribution of lift at subsonic speeds involves basically a 
relatively straightforward process: a system of boimd and 
trailing vortices is set up in the plane of the wing according 
to the specified distribution of lift, and the corresponding 
vortical velocity is calculated, by the Biot-Savart law, at 
points on the surface where the local slopes are desired. 
Keasonably practical numerical and graphical procedures 
have been developed for perfor ming this integration of the 
velocity due to this distribution of vortices (see, for example, 
ref. 1). If the chordwise loading is specified to be uniform, 
as in a number of recent wing-fiesign studies, the problem is 
basically simplified ; as wfil be diown, the solution can then be 
reduced from a double integral over the wing area (or over the 
\ving area plus wake area) to a line integral aroimd the 
boundary of the wing and, in the simplest cases, it can even 
be reduced to a purely analytical procedure. 

The purposes of the present report ore to outlme the basic 
theory behind the solution of these problems involving uni- 
form chordwise loading, to summarize the mathematical 
application of the theory and the development of the required 
formulas, and to describe the actual use of these derived 
results in the design of mean camber surfaces for this type of 
loading. 


The basic theory of the uniformly loaded lifting surface is 
reviewed first. The particular case of the infinitesimally 
wide, uniformly loaded longitudinal strip is next discussed, 
together with the integration of such strips to form the wing 
of arbitrary contour and arbitrary spanwise loading. For the 
unif ormly loaded polygonal wing, closed expressions are 
derived for both the local slope of the mean camber surface 
and the local height of the surface (relative to the leading 
edge). Sections of the mean camber surfaces of four wings 
calculated by these methods are also presented. 

SYMBOLS 

X, y, z streamwise, lateral, and vertical coordinates, 
respectively (see fig. 1) 

a;', 'ifi coordinates of vortex element on wing boimdary 

TJ stream velocity 

w vertical velocity induced by unit vortex (positive 

upward) 

p pressure 

p density 

Gl wing lift coefficient 

Cl wing section lift coefficient 

c chord 

c average chord 

ds vortex element (vector) 

q vector from vortex element to point 

M Mach number 

A aspect ratio 

A sweep angle 

Vortex-segment symbols: 

y^)i (^j 2/a) end points of vortex segment; {yil^yf) 

L=yi CSC a 

M~ —y% CSC a. 
l—^x^~~2xL cos a-f-i* 
m=Va^+2rM cos 
8= V(a;i— a;) ^+yi^ 

i=Va^+2/r 

Subscripts: 

LE leading edge 

TE trailing edge 

I, n, m cases I, n, and HE of appendix 


* Supersedes NAOA TN 2903, “netcrmlnatloa of Mean Camber Sorfoces For Wings Having Uniform Chordwise Loading and Arbitrary Spanwise Loading in Snbsonlo Flow” by 
S. KatiofT, hL Frances Faison, and Hugh C. DuBose, 1963. 
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BASIC THEORY FOR UNIFORMLY LOADED WINGS 


Application of acceleration potential. — flow fields con- 
sisting of a small perturbation flow superimposed on a unifoim 
flow, the pressure is a potential (multiplied by —1/p, it is 
frequently termed acceleration potential) that satisfies 
Laplace^s equation (see, for example, ref. 2, pp.. 225-227). 
In the field of a uniformly loaded lifting surface, then, the 
pressure (relative to free-stream pressure) is a harmonic 
potential that must satisfy the following boundary condi- 
tions: 

(a) It has a uniform negative value over the upper face 
of the lifting surface. 

(b) It has a numerically equal, uniform positive value 
over the lower face of the lifting surface. (That the upper- 
and lower-surface pressures are equal and opposite is not, 
perhaps, obvious merely from the fact that a pressure differ- 
ence exists across the surface. If, however, the lifting sur- 
face is represented by a distribution of boimd and trailing 
vortices, as in ref. 1, this fact is immediately apparent.) 

(c) It vanishes at infinity. 

These boimdary conditions, which uniquely define the pres- 
sure throughout the field, are recognized as identical with 
the conditions on the velocity potential in the field of a 
closed vortex that coincides with the edge, or boimdary, of 
the lifting surface. Accordingly, the pressure at any point 
in the field of a uniformly loaded lifting surface is equal in 
value to the velocity potential of such a vortex, the strength 
of which is the pressure difference between the upper and 
the lower faces, or the lift per unit area. Correspondingly, 
the pressure gradient at any point in the field is equal in 
both magnitude and direction to the potential gradient (that 
is, the velocity) associated with this vortex at that point; 
and it can accordingly be determined by the Biot-Savart 
law. For present purposes, only the vertical component of 
this gradient is of interest. The vortex should not, of 
course, be confused with the lifting vortices of the usual 
airfoil theory; the latter vortices are not used in the present 
report. 

The vertical acceleration of a fluid particle is — ^ ^ (see 

fig. 1 for coordinate system), so that the vertical velocity 
acquired by a particle which has come into the neighborhood 
of the Aving from a large distance upstream is the integral gf 

this expression with respect to time, or 


r® 1 Sp dx 


where the factor dxjU is the element of time. (In^ order to 
simplify the notation, the same symbol x is used for both 
the running variable and the upper limit.) All perturbation 
velocities are assumed te be small so that the path of 
integration, or the path of the fluid particle, is the line 
Constant, Constant. 

Dividing this vertical velocity by Z7 gives the vertical 
slope dzjdx of the streamline. In particular, if the integral 
is evaluated for a point on the lifting surface itself, the local 

1 r* 

slope of the surface is given by — ^ dx, where the 

integration is along the line Constant, 2 = 0 . Here again 
the small-perturbation theory assumes that aU displacements 
from the straight undisturbed streamlines are so small that 
the path of integration may, with sufficient accuracy, be 



Figure 1. — Element of vortex on wing boundary showing coordinate 
and vector systems used in the application of the Biot-Savart law. 


taken in the plane a=0, and, in particular, it assumes that 
the vertical displacement of the trailing edge relative to the 
leading edge is so small that the boundary vortex may also 
be taken in the plane 2 = 0 . The local height of the surface 
2 , relative to the leading edge, is the integral of this slope, or 



For any lift coefficient Ol the pressure difference across 

the lifting surface is pZP, which, as previously noted, is 

numerically equal to the strength of the vortex that is as- 
sumed around the edge of the projected plan form of the 
lifting surface. Thus, finally, the local slope of the surface 
is 

■ ( 1 ) 

where w is the vertical velocity (positive upward) in the 
plane of the lifting surface induced by a unit vortex along 
the edge of the surface. The local height of the surface, 
relative to the leading edge, is then 



The direction of rotation of the unit vortex is such that its 
flow is upward through the surface of the wing; that is, the 
potential increases by unity along a path from the upper to 
the lower wing surface around the edge of the wing. 

Line integral for local slopes of uniformly loaded wing 
with arbitrary plan form. — ^By equation (1), the local slope 
dzjdx may be determined by evaluating the vertical velocity 
w induced by the entire boundary vortex and then integrat- 
ing w from — a> to X. A more convenient method, howevei*, 
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is to evaluate the contribution to w induced by an infinites^ 
mal element of the bounding vortex, to integrate this con- 
tribution from — 00 to a? (which is readily done analytically), 
and then to integrate this result over all the elements of the 
bounding vortex. 

The differential form of the Biot-Savart equation (see ref. 
2, p, 167) for the induced velocity dw due to an element ds 
of a unit vortex is 


dw= 


1 dsXq 

4ir |q|3 


where q is the vector from the vortex element to the point 
in question, and the direction of ds is taken as the direction 
of advance of a right-hand screw rotating in the direction of 
the circulation about ds (see fig. 1). For the present problem, 
where the point hes in the plane of the wing, this induced 
velocity is in the ;^-direction and is given by 


,1 (y—y')dx'~{x—x^dy' 

where (Xf y) is the point at which the induced velocity is 
desired, and (a;', y') is the location of the vortex element 
(da;', dy') on the wing boundary. 

Accordingly, by equation (1), the contribution of a bound- 
ary element (da;', dy') of a uniformly loaded wing to the 
slope of the wing surface at point (a:, y) is 


d 


(Slope) idx*,dv*) 2 ^ 



(2/— ^da/— (x— a/)d2/ 


r (x— a/)da/ . da/ . 

\jy~^y') 

dy' n 

V(x— xO"+(2/— 


The signs of da;' and d'lf are determined from the previously 
mentioned convention for the direction of the vector ds (or 
(da*', d'\f))\ for example, both da/ and dy^ are negative for 
the leading-edge element shown in figure 1. 

The net slope of the mean camber surface at the point 
(Xj y) is the integral of the preceding expression around the 
^ving boundary, or 

xr (a;— a/)da;' 

dx Stt y \jy—y^ 

da/ ^ “I (3) 

where the counterclockwise direction of the integration 
automatically takes care of the signs. The problem of 
determining the local slope of the mean camber sui*face at 
the point (a*, y) is thus reduced to the evaluation of this line 
integral. 


WINGS WITH ARBITRARY PLAN FORM AND 
ARBITRARY SPANWISE LOADING 

Wing considered as sum of uniformly loaded chordwise 
strips. — ^For the wing having uniform chordwise loading and 
arbitrary spanwise loading, it is convenient to consider the 
wing to be made of a series of uniformly loaded chordwise 
strips of infinitesimal span. For each such strip (span dy', 
see fig. 2), the pressure field can be represented by tie 



Figxjbb 2. — Wing composed of uniformly loaded chordwise strips of 
span dy' with a closed vortex superimposed on the boundary of 
each strip. 


velocity potential of a closed vortex superimposed on the 
boundary of the strip. Each of these bounding vortices has 
strength equal to the local pressure difference Ap between 
the upper and lower surfaces of the strip. 

If, as in the preceding analysis, the spanwise loading is 
uniform, aU these closed vortices will be of the same strength, 
so that the chordwise segments common to any two adjacent 
strips cancel and only those vortex elements lying on the 
boundary of the wing remain. The result is thus the same 
as that previously discussed for the uniformly loaded wing 
(eq. (3)). 

If, however, the spanwise loading is not uniform, the 
closed vortices surrounding adjacent strips will be of tmequal 
strength; vortex segments common to adjacent strips will 
no longer cancel, and vortex elements lying on the wing 
boimdary will vary in strength along the boimdary. The 
contribution of these boundary vortex elements to the in- 
duced velocity w can stUl be summed by a line integration 
around the wing boimdary of the expression given in equa- 
tion (3), except that Cl must be replaced by the local section 
lift coefficient Ci(y^) and placed under the integral sign. The 
contribution of the chordwise segments is derived in the 
following paragraph. 

Streamwise vortex segments. — ^The Biot-Savart formula 
for the induced velocity due to the straight-line chordwise 
vortex segment of unit sti*ength is 

14;==—^ (cos 02 +cos ^i) 


( 4 ) 
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where r, di, and 6$ are defined in figiure 3. With r and the 
cosines expressed in Cartesian coordinates, this expression 
becomes 

1 ( a;— (a/+c) z—x' ) 

^ 4v(y—y^\-\l[x—ix'+c)Y+(y—yO^ -^{x—xTf+iy—yy] 

where c is the local diord. The integral of this expression 

from — ® to ar, multiphed by — dcj/2 or — ^ ^ dy' (see eq. 

(1)), gives the contribution to the slope from the streamwise 
vortex segment of strength corresponding to dci and length 
equal to the local chord: 

d(Slope)= — (y— ^ ^ (tVC®— ®0*+(2/— yO*+(a:— — 

y{x-ix'+c)V+(y-yO^+x-ix'+c)]) 

Finally, integrating tliis last equation with respect to y' 
across the span of the wing (from left to right) gives the 
contribution of all these chordwise segments to the slope at 
the point It will be observed, however, that the term 

within the second bracket in this equation is the same as 



Figure 3. — Geometrical relationships for straight-line chordwise vortex 
segments that occur with spanwise-varying area loading. 


that within the first bracket except that a;' is replaced by 
x'+c, the corresponding tr ailing-edge coordinate. Accord- 
ingly, the integral of this expression across the span can, for 
convenience, be considered as the line integral 

[V(a;-®0*+(2/-2/0*+(a:-a:')] (5a) 

or 

^ ^ [V(*— 0"+(y-2/')"+(a:— aO] (5b) 

around the boundary of the wing. 

The problem of determining the local slope of the mean 
camber surface at point (a;, 2 /) is thus reduced to the evaluation 
of the line integrals (3) and (5), where the integral (6) is 
omitted if there is no spanwise variation in the area loading 
(or in the local lift coej0&cient Cj) and where the local lift 
coefl&cient Ci replaces Oi, in equation (3) and must be brought 
imder the integral sign if there is a spanwise variation of Cu 

COMPUTATION 

Although computing the integrands of expressions (3) and 
(5) and then performing the integrations should be a fairly 
straightforward process, a short outline of suggested proce- 
dures may be helpful. It may be noted at the beginning 
that, since the slopes (and the integrals for the slopes) are 
nondimensional, the results will be independent of the dimen- 
sional scale chosen for the work; taking the root chord or 
the semispan as unity will probably be most convenient. 
It may also be remarked that a carefully drawn plan form of 
the wing wiQ be helpful in setting up the computations. 

ESsumS of procedure for computing slopes of uniformly 
loaded mean camber surface. — Given the plan form of a 
wing that is to have uniform area loading and a specified lift 
coefl&cient, a possible procedure is as follows: 

(1) Select points (jjy) on the surface where the slopes are 
to be obtained. In general, these points should lie along 
several selected chord lines, with perhaps four along each 
line. 

(2) Select points (x',y') along the leading and trailing 
edges where the integrands of equation (3) are to be deter- 
mined. In gener^, trading-edge points should bo at the 
same spanwise positions as the leading-edge points. 

(3) Consider the integral in equation (3) to be broken up 
into two parts — one with respect to x' and one with respect 

and 

!^V(*-aO*+(2/-J/0^ 

For each point (a:, 3 /), compute the values of the integrands in 
expressions (6a) and (6b) for all the (x',y') points, 

(4) For each point lx,y), plot the integrands of (6a) against 
x', plot the iutegrands of (6b) against y', and determine the 
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area of each (for example, by ru nn ing the planimeter aroxmd 
the curve in the direction corresponding to a counterclock- 
wise movement of the variable point around the wing 

boundary). 

(5) For each point {x,y), sum these two areas and multiply 
by — Cx/Stt in order to determine the local slope dzjdx, 

(6) The actual heights of the mean camber surface along 
the chosen chords are now determined by integrating the 
slopes foimd in the preceding steps. The integration cannot 
extend quite to the leading and trailing edges, however, 
because the slopes cannot be readily determined very close 
to the ends of the chord lines, where the slopes become 
infinite. Near the ends, however, the mean camber lines of 
the sections may be considered as NACA a=1.0 mean lines 
(ref. 3), and they may accordingly be filled in, with generally 
adequate accuracy, by fitting NACA a— 1,0 mean lines 
between the front and rear limits of the calculated segments 
of the mean camber lines and the leading and trailing edges. 

The integral (6a) is actually improper, since the integrand 
becomes infinite where y'=y. The Cauchy principal value 
exists, however, and is evaluated by first adding the inte- 
grands at equal distances but on opposite sides from the 
singulai* point and then integrating the sum. This method 
will be further described in a subsequent section. 

R6sum6 of procedure for computing slopes of a mean 
camber surface with uniform chordwise loading and arbi- 
trary spanwise loading. — ^For an arbitrary spanwise loading, 
an additional integration is required, together with a modi- 
fication of the preceding integrations. 

(1) Determine the integrals 

and 

by the same process as was used for (6a) and (6b), except 
that Cif being now a variable, must be brought inside the 
integral. 

(2) Determine the spanwise rate of change of local lift 
coefficient dcijdy' at points y\ Compute values of the 
integrand in 

Plot against y' and integrate. This integral is also improper 
and is treated as previously mentioned. 

(3) For each point (a?,i/), the local slope dzjdx is 

^ [Integral (8)— Integral (7a)— Integral (7b)] 

07T 

Example. — In figure 4 is shown the plan form of a swept 
^ving and the desired spanwise lift distribution. As an 


example of the computation, integrals (7a), (7b), and (8) 
will be obtained for the point designated P in the figure. 
Points designated a . . . z and a' . . . z' are the points on 
the boundary where the integrands were evaluated; the 
primed symbols are used merely because more than 26 sym- 
bols were needed. The origin was taken at the wing apex 
and the wing semispan was assumed equal to unity. 

The computation of the integrands at the points around the 
boundary is given in table I. Most of the points were evenly 
spaced, but near the singularity (y'=y), the intervals were 
reduced to one-tenth as much as most of the other intervals. 

Figure 5 shows the plots of the integrands (multiplied by 
1/87T or — 1/87 t) against a;' or y^ and also shows the method of 
determining the Cauchy principal value at the singularity. 
For example, in figure 5(c) it will be seen that the integrand 
goes to 00 at the right of the singularity and to — oo at the 
left of the singularity. Algebraically adding the integrands 
at equal distances from the singularity (for example, the 
value at point b plus the value at point v, the value at point 
c plus the value at point u, and so on) and plotting the sum 
results in the section B1 of the curve, where the portion 
nesirest the singularity is obtained by extrapolation from 
point k. The desired integral (8) is finally determined by 
nirming the planimeter along the path abBILVva'b'jVs'a. 

As a further example of the intermediate steps in the 
calculation of a mean camber surface, figure 6(a) is shown, 
which is a plot of the slope dzjdx along the streamwise chord 
throigh point P. The example calculation described in the 
preceding paragraphs, it will be noted, gives the three in- 
tegrals the sum of which provides one point on this dzjdx 
curve of figure 6(a). In figure 6(b) is shown the correspond- 
ing curve for the NACA a=1.0 mean line, for C|=1.0, from 
reference 3. The curve may be found useful in extra- 
polating to the leading and trailing edges, as previously 
mentioned. 

As may be inferred from the preceding rfeum6 of the 
computations for one point P, the total effort required to 
compute accurately a mean camber surface by the method 
described is very large, although, according to the authors’ 
experience, it is not at all prohibitive. Nevertheless, the 
work is ideally suited to modem high-speed computing 
machinery — ^for example, of the punched-card type — so that 
consideration should be given to the use of such equipment 
where it is available. 

It may also be mentioned that two contour charts of the 
integrands in (6a), (6b), and (8), with the factor Ci or dcijdy 
omitted, have been prepared on film and may be obtained on 
request from the National Advisory Committee for Aero- 
nautics. Their form is such that, if the transparency is 
superimposed on a correctly scaled drawing of the wing 
plan form, with a boundary point {x\y') at the origin of the 
chart, the contour value read at point (x,y) is the desired 
value of the integrand. The charts are satisfactorily 
accurate except in the neighborhood of the singularities. 
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Figure 4. — Plan form and spanwise lift distribu 
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TABLE I 

COMPUTATION OF INTEGRANDS FOR INTEGRALS (7a), (7b), AND (8) 
FOR POINT P (^=0.7707, y=0.2600) OF FIGURE 4 


Boundary point, 

x—x* 

<D 

y-yf 


d 

dcijdy' 

1 — 1 

® et© 

1 I 

®|0 


(?)••• 


a 

0. 7707 

0. 2600 

0. 8102 

'0. 6880 

0 

4. 5900 

0. 7259 

0 

b 

.4660 

. 1260 

.4826 

. 6725 

-.214 

9. 0097 

1. 1874 

-1. 6238 

c 

.4355 

. 1125 

.4498 

. 5700 

228 

9. 9722 

1. 2672 

-1. 7942 

d 

.4051 

. 1000 

.4173 

. 5669 

-.244 

11. 1726 

1. 3586 

-2. 0087 

e 

. 3746 

.0875 

. 3847 

. 5637 

-.259 

12. 7161 

L4663 

-2. 2475 

f 

.3441 

.0760 

. 3522 

. 6604 

-. 273 

14. 7721 

1. 5911 

—2. 6345 

g 

. 3137 

.0626 

. 3199 

. 6669 

289 

17. 6479 

L7409 

—2. 9298 

h 

. 2832 

. 0500 

. 2876 

. 5536 

300 

21. 9746 

1. 9249 

-a 4248 

i 

.2527 

. 0376 

. 2555 

. 6601 

-. 316 

29. 1773 

2. 1530 

-4.2824 

j 

.2223 

.0250 

.2237 

.5459 

-. 331 

4a 5344 

2.4403 

— a 9050 

k 

. 1918 

. 0126 

. 1922 

. 6419 

-. 344 

8a 6130 

2. 8195 

-10. 6677 

1 

. 1613 

0 

. 1613 

. 5373 

-. 360 


a 3329 


m 

. 1309 

-. 0126 

. 1316 

. 6329 

-. 373 

-85. 0679 

4.0525 

7. 8300 

n 

. H)04 

-. 0250 

. 1035 

. 5279 

-. 387 

-41. 6986 

a 1006 

a 1564 

0 

.0699 

-. 0376 

. 0793 

. 6226 

-.401 

-26.2156 

a 5889 

1. 5954 

P 

.0394 

-. 0500 

.0637 

. 6168 

-.415 

-16. 7288 

a 1130 

.8557 

q 

. 0090 

-. 0625 

.0631 

. 5115 

-.429 

-9. 3610 

a 1062 

.4949 

r 

0215 

-. 0750 

.0780 

. 5059 

-.444 

-4. 8863 

a 4869 

. 3345 

s 

-. 0520 

0876 

. 1018 

. 6009 

-.458 

-2. 8005 

4. 9204 

.2607 

t 

-. 0824 

-. 1000 

. 1296 

.4947 

-. 470 

-1. 8017 

a 8171 

.2218 

u 

1129 

-. 1125 

. 1694 

.4885 

-. 484 

-1. 2666 

a 0646 

.2001 

V 

-. 1434 

-. 1260 

. 1902 

.4831 

-. 498 

-. 9497 

2. 5363 

. 1865 

w 

-.4481 

-. 2500 

. 6131 

. 4123 

-. 634 

-. 2091 

.8039 

. 1648 

X 

-. 7527 

-. 3760 

.8409 

. 326L 

-. 770 

-. 0912 

.3879 

. 1811 

y 

-1. 0574 

-. 5000 

1. 1697 

. 2206 

-. 902 

— . 0426 

. 1890 

.2026 

z 

-1. 3621 

-. 6260 

1. 4986 

. 1006 , 

-1. 006 

-. 0146 

.0668 

.2195 

a' 

-1. 6668 

-. 7500 

L 8278 

0 

-. 380 

0 

0 

.0816 

b' 

-2. 5271 

7500 

2. 6360 

0 

-. 380 

0 

0 

.0552 

c' 

-2. 2941 

6260 

2. 3777 

. 1006 

-1. 006 

-. 0066 

.0421 

. 1344 

d' 

-2. 0611 

6000 

2. 1209 

. 2206 

-. 902 

-. 0125 

. 1042 

. 1079 

e' 

-1. 8281 

-. 3760 

1. 8662 

. 3261 

-. 770 

-. 0177 

. 1748 

. 0782 

f' 

-1. 6951 

-. 2500 

1. 6146 

.4123 

-. 634 

-. 0200 

.2555 

. 0495 

g' 

-1. 3621 

-. 1250 

1. 3678 

.4831 

-.498 

-. 0162 

. 3527 

. 0227 

h' 

-1. 1291 

0 

1. 1291 

. 5373 

-. 360 

0 

.4761 

0 

I' 

-. 8961 

. 1260 

. 9048 

. 6726 

-.214 

.0439 

. 6322 

-. 0149 

j' 

1 -. 6631 

. 2500 

. 7078 

. 6880 

0 

. 1512 

.8297 

0 

k' 

-. 8961 

. 3750 

. 9714 

. 5725 

.214 

. 1182 

. 5888 

, 0432 

V 

1 -1.1291 

. 6000 

1. 2349 

. 6373 

. 360 

. 0920 

.4348 

. 0760 

m' 

-1. 3621 

. 6250 

1. 4986 

.4831 

.498 

.0705 

. 3228 

. 1088 

n' 

-1. 6951 

. 7500 

L 7626 

.4123 

. 634 

. 0622 

.2337 

. 1414 

o' 

-1. 8281 

.8760 

2. 0267 

. 3261 

.770 

. 0366 

. 1608 

. 1748 

P' 

-2. 0611 

1. 0000 

2. 2909 

. 2206 

.902 

.0221 

1 . 0960 

.2070 

q' 

-2. 2941 

1. 1250 

2. 6651 

. 1006 

1. 005 

. 0092 

. 0395 

. 2329 

r' 

-2. 5271 

L2500 

2. 8193 

0 

. 380 

0 ■ 

0 

. 0888 

s' 

-1. 6668 

1. 2600 

2.0834 

0 

. 380 

0 

0 

. 1266 

t' 

-1. 3621 

1. 1250 

1. 7666 

. 1006 

1. 005' 

.0206 

. 0572 

.3610 

u' 

-1. 0674 

1. 0000 

1. 4564 

.2206 

. 902 

. 0602 

. 1512 

.3586 

v' 

-. 7627 

.8750 

1. 1542 

. 3261 

. 770 

. 1296 

. 2824 

. 3533 

w' 

-.4481 

. 7600 

.8737 

. 4123 

. 634 

.2676 

.4716 

. 3600 

x' 

-. 1434 

. 6260 

. 6412 

.4831 

.498 

. . 6010 

. 7646 

.3966 

y' 

. 1613 

. 6000 

. 5254 

. 6373 

. 360 

1. 4035 

1. 0219 

.4931 

z' 

. 4660 

. 3750 

. 5981 

. 6725 

.214 

.2. 7137 

. 9564 

. 6101 


♦Integrand for (7a) 
♦♦Integrand for (7b) 
♦♦♦Integrand for (8) 




SrO'-y') ^(x-xY + (y-y')^ 
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(a) Calculated slopes, example wing. 

FiauBB 0. — Calculated slopes along streamwise mean camber line 
through point P (see fig. 4) and slopes along NACA a=1.0 mean 
lino at c/=1.0. 


POLYGONAL WINGS WITH UNIFORM AREA LOADING 


For a uniformly loaded wing (uniform chordwise loading 
and uniform Ci across the span) the boundary of which con- 
sists of a munber of straight-line segments, the integrations 
to determine the local slope of the surface and the local 
height of the surface may both he done analytically. The 
previously derived equations can be integrated with respect 
to a?' and y', where the point (a?',y') moves from one end of 
the straight-line segment to the other. In the present 
development, however, the straight-line segment will be 
treated as a whole. The vertical induced velocity w due to 
the unit bounding vortex is expressed as the sum of Wa.^ 
v)bj . . the velocities due to the separate straight seg- 
ments A, B, . , . (see fig. 7) which are given by the Biot- 
Savart law (eq. (4)). The contribution of segment A to the 


, slope of the surface at point {x,y) is, therefore. 



and the contribution to the height of the sxnface at the 
point (a?,y) is 



w^dx 



(b) Slopes along NACA a=1.0 mean line, for c/=1.0, from reference 3. 
Figuhb 6. — Concluded. 


Summing these expressions for aU the segments A, B, 

. gives the total slope or height of the mean camber 
surface at the desired point. 

As already noted, these integrals for the separate straight 
segments can be evaluated analytically. Because the 
mathematical manipulation and the resulting formulas are 
somewhat lengthy, they are given in the appendix. Tliree 
different cases, distinguished by the relative geometry of 
the vortex segment (that is, the segment of the wing bound- 
ary) and the point {x^y) where the slope or height of the 
surface is to be found, are discussed in the appendix. In 
case I (fig. 8), the path of integration from — to aj crosses 
the segment; in case II (fig. 9), the path of integration does 
not cross the segment; and in case HI (fig. 10), which is a 
special case of II, the path of integration is parallel to the 
segment. 

It should be noted that this same problem has been treated 
from a somewhat different viewpoint in reference 4. 

COMPRESSIBILITY CORRECTION 

If the mean camber surface is desired for a compressible 
subsonic flow at Mach number Af, the Prandtl-Glauert 
method, as described in reference 5, may be used. That is, 
aU the longitudinal (streamwise) dimensions of the wing are 
stretched by the factor l/Vl’-ibP, so that the aspect ratio is 
reduced by the factor and the tangent of the sweep 

angle is increased by the factor 1/Vl—ibf^ and the mean 
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Figure 7. — Unit bounding vortex composed of separate straigtit 
segments A, B, C, D, E, and F superimposed on vring plan form 
boundary. 

camber surface is i^alculated for this fictitious wing in in- 
compressible fiow at a lift coefficient equal to the desired lift 
coefficient multiphed by 1— ikP. The ordinates (values of 2 ) 
so obtained wiQ be precisdy the ordinates of the mean 
camber surface for the desired wing at the desired lift cp- 
efficient in the compressible fiow at Mach number M, 

EXAMPLES OF CALCULATED MEAN CAMBER SURFACES 

Sweptback wing, .4=1,74. — ^Figure 11 shows (by the 
curved lines) mean camber lines calculated for the wing 
already described in the example and in figure 4. The points 
where the slopes were computed are indicated by small 
circles. It should be noted that figure 11 and subsequent 
similar figures do not represent oblique projections of the 
wings. Rather, they show the true plan forms of the wings 
(in the a^-plane) and the true mean camber lines (parallel to 
the aJ 2 -plane) ; accordingly, the y- and 2 -axes coincide on the 
figures. 

In figure 11 and the subsequent similar figures, aU the 
leading-edge points have been assumed to lie along a hori- 
zontal line. This choice is, of course, arbitrary; that is, the 
extent to which the present linear theory is applicable would 
be essentially unaltered if, for example, the wing had a 
reasonable amount of dihedral. 

Uniformly loaded triangular wing, — The formulas derived 
in the appendix for cases I and II were used to calculate the 
mean camber surface of a triangular wing having an angle 


of sweep of 68.4° of the leading edge (aspect ratio, 1.67) such 
that the wing should be uniformly loaded at unit lift co- 
efficient in incompressible flow. Figure 12 shows the plan 
form and the chordwise camber lines for several spamviso 
stations. The dimensions shown correspond to a root chord 
of unity; however, the 2 -scale is shown somewhat exaggerated. 

Uniformly loaded swept wing. — The formulas derived for 
cases I, n, and.m were used to calculate the mean camber 
surface of a swept wing of aspect ratio 8, taper ratio 0.46 
and 45° sweepback of the quarter-chord line such that the 
wing should be uniformly loaded at unit lift coefficient at a 
Mach number of 0.9. In accordance with the proposed method 
of taking into account compressibility, the calculations were 
made for the stretched wing in incompressible flow. The 

stretching factor is . __ — =2.29, so that the aspect ratio 

Vl-lO.97 

of the stretched wing is 3.5 and the sweep of the quarter- 
chord fine is 66.46°. Figure 13 (a) shows this stretched ^ving 
and the mean camber lines calculated for this stretched Aving 
in incompressible flow at unit lift coefficient. The dimen- 
sions shown correspond to a semispan of unity; however, 
the 2 -scale is shown somewhat reduced. The corresponding 
mean camber lines 'for the physical Aving (ri=8; design 
(7i=1.0 at M=0.9) should have 44 percent as much camber 
as the mean camber lines for this stretched Aving in the incom- 
pressible flow. 

Swept vdng with eUiptioal span load distribution. — For the 
wing plan form of the preceding example another mean 
camber surface was calculated that, at unit lift coefficient and 
a Mach number of 0.9, would give uniform chordAvise loading 
but an elliptical span load distribution for the Aving as a 
whole. Mean camber lines for this case are presented in 
figure 13 (b). 

This example was also calculated, under the direction of 
Mr. Robert R. Graham of the Langley Laboratory, by the 
method of reference 1. The two results were in very good 
agreement. No definite information was obtained, hoAvever, 
Avith regard to the relative expediency of the tAvo methods 
(that of ref. 1 and that of the present report). One might 
suppose that the present method would be preferable for 
Avings Avith uniform chordAvise loading, since it is designed to 
take advantage of this particular characteristic. The Avork 
represented by table I and figure 5 (outlining the computa- 
tions for the slope at one point), however, is by no means 
small, so that such a presumption is not de^tely sub- 
stantiated by present experience. Perhaps the fact that 
the computations and integrations are of such form that 
they can be readily performed by modern high-speed com- 
puting machinery constitutes the most significant charac- 
teristic of the present method. 

Langley Aeronautical Laboratory, 

National Adausory Committee for Aeronautics, 
Langley Field, Va., January IS, 195S. 
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Figure 8. — Geometrical relationships for case I. Heavy lines in 
small sketch show relative positions of vortex segment and path of 
integration for this case. 



Figure 9. — Geometrical relationships for case II. Heavy lines in 
small sketches show relative positions of vortex and path of integra- 
tion for three different conditions for which this case applies. 
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(a) Uniform area loading. 

(b) Uniform chordwise loading; elliptical span load distribution. 

Figueb 13. Mean camber lines for the stretched wing at unit lift coeflGlcient in incompressible flo\V. Aspect ratio, 3.6; taper ratio, 0.45; quarter- 
chord sweep angle, 66.46®; Ol— 1.0. For the physical wing (aspect ratio, 8; taper ratio, 0.46; quarter-chord sweep angle, 45®) with unit lift 
coefficient at il/=0.9, the mean camber lines should have 44 percent as much camber as the lines shown here. Small figures show span load 
distributions. Circles indicate points for which z was computed. 





APPENDIX 

DEVELOPMENT OF FORMULAS FOR UNIFORMLY LOADED POLYGONAL WINGS 


Case I — Patli of integration crosses vortex segment. — If 
the vortex segment is the leading edge of the ving, the path 
of integration may cross it. The vortex segment and the 
path of integration for this case are shomi by the heavy 

lines in the small sketch in figure 8. The integral j v)idx 

is improper because the integrand becomes infinite where 
the path of integration crosses the vortex. The Cauchy 
principal value, however, can be determined. Before the 
integration is performed, the Biot-Savart formula (eq. (4)) 
is expressed in terms of the variable x and the fixed dimen- 
sions of the vortex segment, where the origin is defined as 
the point of intersection of the path of integration and the 
vortex segment. Accordingly, the end points of the vortex 
segment are (a?i,2/i) and ^, 2 / 2 ), where t/i> 2/2 (see fig. 8). Let 

X=2/i CSC a 

7/2CSC a 

Then it can be seen that 


cos ^ 1 = 


r—x sm a 

Z—x cos a L—'X cos a 


cos 02 = 


^ 2xL cos a+Z® 

M-\-X cos oc M-\-X cos CL 


m 


Va^+2xM cos a+i\P 


The contribution of Wi to the slope of the mean camber 
surface will be 


=_i 0 r 1 i 


Z — X cos a 


2xL cos a+L^ 


M-{-x cos a 


cos 

-o^ 


Sir sm 


log.[j 




+ 


x*sin* 


Q.+L—X cos a) (m+M+x cos 


«)]} 


(Al) 


It is of interest to note that, along the wing tip or the wing 
root, where L or M goes to zero while x is positive, this slope 
becomes infinite. 

The contribution to the height of the mean camber surface 
wiU be the integral of this last expression; that is, 
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dx 




{ J2 

X cos a log, ]jj2+2x log, sin a+x log, x®+ 


Stt sin I 

(x cos a—M sin^ a) log, (x+M cos a+^) — 

(x cos a+L sin^ a) log, (x— Z cos — 

X log, {m+M+x cos a)(Z-t-Z— X cos a)+cos a(Z— m)^ | 

(A2) 

or, after substituting the indicated limits, 

Q r 

gi= 8^ sin a I ® log, sin a+x log, x^-4- 

(x cos a—M sin^ a) log, (x+M cos a+m) — 

(x cos a+Z sin^ a) log, (x— Z cos a+Z) — 

X log, (m+M+x cos a)(Z+Z— X cos a) + 

M sin^ a log, M(l+cos «)+ 

Z sin* a log, Z(l— cos a)+I(Z— m)— (Z— M)] cos 

(AS) 

Case n — Path of mtegration does not cross vortex seg- 
ment. — ^The three small sketches in figure 9 show thi’eo cases 
in which the path of integration does not cross the vortex 
segment. In two of these cases, the segment hes along the 
trailing edge; in the third case, the segment lies along the 
leading edge but hes wholly to one side of the path of inte- 
gration. 

For the derivation of the formulas for case II, the origin of 
coordinates is specified , as the intersection of the path of 
integration with the vortex or with the line of the vortex 
extended. The end points of the vortex segment are again 
designated (xi,yi) and (xa,^,), where yi>^ 2 . The angle a is 
defined as before. Following the sign convention indicated 
in figure 8, M will be considered positive if (x 2 , 2 / 2 ) hes to the 
left of the origin, and Z will be considered positive if (xi,yi) 
hes to the right of the origin. For example, in figure 9, Z is a 
negative quantity, whereas M is a positive quantity. 

For any of the three cases shown in figure 9, the expression 
for 'Wn is identical with that for Fm-thermoro, the ex- 
pression for 




is identical in form ■with the expression (Al) derived for 

j that is, neither singular points in the integration nor 

negative values of x, L, or M affect the expression for the 
integral. 
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Intograting ^ order to obtain Zn yields the same 

expression for the indefinite integral as was found in inte- 
(I), ; SO that the only diJfferences between the ex- 
pressions for Zii and Zi arise from the different hmits of iute- 
gration— for the limits were 0 and x, whereas for Zu, the 
limits are Xlb a^nd x; that is, 


Zji' 


X'XS)./" 

Q r j2 

=s -< a; cos a log, ■^+2a: log^ sin a-\-x log. a^+ 

Sir sm a JM 

{x COS a—M sin^ a) log« (x+M cos a+m) — 

(x cos a+i sin^ a) log, (x— i cos a+Z) — 

xlog, (m+M+a; cos a)(Z+L—a;cos a)+cos o:(Z— m)^j^ 


rLE 

(A4) 


Case m — ^Vortex segment parallel to free stream and 
hence to path of integration (fig. 10). — Coordinates are 
chosen so that the path of integration lies along y=0. The 
end points of the vortex segment are (0, yO and (xi, where 

a:t>0. 

Then, from figure 10, 

s=-^(xi—xY+yi‘ 

t=V^+^ 


cos 01 = 


Xi—X Xi — X 

8 -yjixi—xy+yi^ 


. X X 

cos 0 a= 7 = 


^ Va^+2/i* 


Then 

^ r i~ I ^ 

Siryi J_„ l_Y(a;,— a;)='+2/i* -^x?+yi^ 

=— yx^+yi^—^{xi—xy+yi~+xi] (A5) 

and 

2ra=— ({(a;i— x) ^|iXl—xy+yy+x■^J3^+yy+ 

yyiog, [xi — x+ V (a:i— x)*+i/i“J + 

yy log, (x+ V x*+ 2 /t*)}— foi— aifi) -J{xi—XLE)-+yy+ 

Xle ^/xLB^+yy+yyiog, [xi— XiB+Vfe— 

yy log. [xi«+Vais*+3/i*+2xi(x— Xis)!}) (A6) 

For equations (A5) and (A6), the direction of rotation of 
the vortex segment was assumed to be that corresponding 
to the right wing tip. For the left wing tip, the signs 
should be reversed. Stated differently, the equations will 
be correct iu either case if 2/1 in the first factor on the right- 
hand side is replace by I 2 / 1 I. 
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